Itinérant Ferromagnetism in ultracold Fermi gases 
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Itinérant ferromagnetism in cold Fermi gases with répulsive interactions is studied applying the 
Jastrow-Slater approximation generalized to finite polarization and température. For two compo- 
nents at zéro température a second order transition is found at akF — 0.90 compatible with QMC. 
Thermodynamic functions and observables such as the compressibility and spin susceptibility and 
the resulting fluctuations in number and spin are calculated. For trapped gases the resulting cloud 
radii and kinetic énergies are calculated and compared to récent experiments. Spin polarized Sys- 
tems are recommended for effective séparation of large ferromagnetic domains. Collective modes 
are predicted and tri-critical points are calculated for mufti- component Systems. 

PACS numbers: 71.10.Ca, 03.75.Ss, 32.80.Pj 
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I. INTRODUCTION 

Ultracold Fermi Systems with strong attraction be- 
tween atoms has led to important discoveries as universal 
physics and the BCS-BEC crossover. Recently strongly 
répulsive interactions has been studied and a transition 
to a ferromagnetic phase was observed in the experiments 
of Jo et al. [H . Earlier Bourdel et al. and Gupta et 
al. H also observed a transition when the interactions 
became strongly répulsive near Feshbach résonances. A 
phase transition from a paramagnetic (PM) to ferromag- 
netic (FM) phase was predicted long ago by Stoner [4| 
based on the Hartree-Foch mean field energy and has 
recently been confirmed by more elaborate calculations 
including fluctuations @,@| and by QMC The cal- 

culated transition points and order of the transition differ 
also from experiment The FM transition is disputed 
by Zhai [9] who claims that the expérimental data is com- 
patible with strongly correlated répulsive Fermi Systems 
which would explain the inability to observe FM domains 
in Réf. fl. 

It the purpose of this work to clarify the phase dia- 
gram of strongly répulsive Fermi atomic Systems as well 
as to calculate thermodynamic functions and measurable 
observables in atomic traps that clearly can distinguish 
the FM and PM phases and détermine the order of the 
transition and the universal functions. By extending the 
Jastrow-Slater model [ÏÔI - [Ï^ | to finite polarization and 
température, we calculate the free energy and find a sec- 
ond order FM transition in a répulsive Fermi gas. A num- 
ber of thermodynamic functions as the spin susceptibil- 
ity, compressibility, and observables as radii and kinetic 
énergies can be compared to experiments, and others as 
fluctuations, collective oscillations and phase séparation 
can be predicted. 

As a start the dilute limit model of Stoner is extended 
to finite température and the polarization and order of 
the transition is determined and compared to second or- 
der calculations. Subsequently, the Jastrow-Slater ap- 
proximation is described for the correlated manybody 
wave-function in the strongly interacting limit and ex- 
tended to finite polarization and température. Detailed 



calculations of the free energy and a number of ther- 
modynamic functions are given. In particular the spin- 
susceptibility and compressibility are used for calculating 
fluctuations in spin and total particle number in section 
III. In section IV finite traps are considered and the cloud 
radii and kinetic énergies are calculated and compared 
to récent experiments Collective modes are discussed 
in section V. Multi-components Systems are discussed in 
section VI and a new string of critical points is found and 
plotted in a multi-component phase diagram. Finally, a 
summary and outlook is given. 



II. FERROMAGNETIC TRANSITION 

The models for répulsive ultracold Fermi gases in Refs. 
ail predict a phase transition somewhere near the 
unitarity limit akp ~ 1 but the phase diagrams disagree 
quantitatively as well as qualitatively concerning the or- 
der and critical points. 

For a référence model we start with a simple finite 
température extension of the Hartree-Fock approxima- 
tion originally studied by Stoner [4], which is a dilute 
limit expansion to first order in the scattering length. 
Subsequently, we calculate the phase diagram in the JS 
approximation and compare to those in the dilute limit 
to first [i| and second [&|(| order as well as QMC 0,9- 



A. Dilute approximations 



A dilute (kpa <C 1) degenerate Fermi gas with atoms 
in spin states a = 1, ..,v with densities p a = k Fa ./6Tr 2 
and Fermi energy Ep^ = ksTp^ = h 2 k F a /2m has the 
free energy 



E4na 
PaPa 
m 



(1) 



It consists of the kinetic energy, the interaction en- 
ergy to lowest order in the scattering length a as 
in the Stoner model [4], and the thermal energy 
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Figure 1: Phase diagram for a two-component Fermi gas with 
répulsive interactions. Full (dashed) curves indicate first (sec- 
ond) order PM to FM transitions within JS and dilute ap- 
proximations with @ and without Q fluctuations. The circle 
indicates a tri-critical point. Triangles show the QMC transi- 
tion points at zéro température of Refs. @, @|- 



f T = —(-ïï 2 /4)(m,*/m)Y^ r jPaT 2 /E F}< j at low tempéra- 
tures r « È^j. In the dilute limit the effective mass 
m* /m = 1 + [8(71n2 — 1) / I5ir 2 ]a 2 k 2 F in two-component 
symmetric Systems only déviâtes from the bare mass to 
second order in the interaction parameter ak F - 

The density of the components are equal only in the 
PM phase and when the components are balanced ini- 
tially. In the following we define an average Fermi wave 
number kp from the total density p = i/kp/Gn 2 . 

We postpone multicomponent Systems to sec. VI and 
concentrate first on two spin states, e.g. a =f, J, with 
total density p = Pi + p^. The population of spin states 
are allowed to change (polarize) in order to observe phase 
transitions to itinérant ferromagnetism. The polarization 
(or magnetization) P — (p^ — Pî)/p of the ground state 
phase is found by minimizing the free energy at zéro mag- 
netic field. The free energy of a low température idéal 
gas is f = E F p(3/5 - i: 2 T 2 /AE 2 F ). Expanding Eq. (TJ 
for small polarization leads to a Ginzburg-Landau type 
équation for the free energy 
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to leading orders in interaction, polarization and tempér- 
ature. Here xt = (d 2 f /dP 2 )' 1 is the isothermal spin- 
susceptibility given by 



X° 1 ^ tt 2 T 2 
— = 1 akp H <r 

XT 12 T 2 



(3) 



where xo — 3p/2Ep is the spin-susceptibility for an idéal 
gas at zéro température, xt becomes singular when 
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where the free energy of Eq. @ predicts a second order 
phase transition from a PM to a FM (see Fig. 1) in 
accordance with the zéro temp érature rcsult of Stoner 
U • The polarization is P = ±y/27{ak F /n - 1/2) at zéro 
température but quickly leads to a locally fully polarized 
System P = ±1 due to the small fourth order coefficient 
inEq. ©. 

However, the predicted transition occurs close to the 
unitarity limit where the dilute équation of state Eq. ((T|) 
is not valid. Higher orders are important as exemplified 
by including fluctuations, i.e. the next order a 2 correc- 
tion. As found in Refs. [!, Q fluctuations change the 
transition from second to first order at low températures 
up to a tri-critical point at température ~ 0.2Tf, where 
the transition becomes second order again (see Fig. 1). 
However, the 2nd order expansion is not valid either in 
the unitarity limit. 



B. Jastrow-Slater approximation 

The JS approximation applies to both strongly attrac- 
tive and répulsive crossovers where it already has proven 
to be quite accurate for predicting universal functions 
and parameters. The JS approximation is the lowest 
order in a constrained variational (LOCV) approach to 
calculate the ground state énergies of strongly correlated 
Systems. It was developed for strongly interacting and 
correlated Bose and Fermi fluids respectively such as 4 He, 
3 He and nuclear matter [lOj. JS was among the earliest 
models applied to the unitarity limit and crossover of 
ultracold Fermi (TU and Bose [Ï2| atomic gases. As ex- 



plained in (ÏÔl4l2l | the JS wave function 



(5) 



incorporâtes essential two-body corrélations in the Jas- 
trow function </>(r). The antisymmetric Slater wave func- 
tion $5 for free fermions $5 insures that same spins 
are spatially anti-symmetric. The Jastrow wave function 
only applies to particles with différent spins (indicated by 
the primes). The pair corrélation function <p(r) can be 
determined variationally by minimizing the expectation 
value of the energy, E/N = {^\H\^) / Wy), which may 



be calculated by Monte Carlo methods [8|, [13j- At dis- 
tances shorter than the interparticle spacing two-body 
clusters dominate and the Jastrow wave function cj>(r) 
obeys the Schrôdinger équation for a pair of particles of 
différent spins interacting through a potential U(r) 



m dr 2 



U(r) 



-) = 2\r<j>{r) , 



(6) 



where the eigenvalue is the interaction energy of one atom 
A = E int /N. Most importantly, the boundary condition 
at short distances (r = 0) is given by the scattering length 



{r4>)> 



(7) 
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Many-body effects become important when r is compa- 
rable to the interparticle distance ~ fc^ 1 , but are found 
to be small |10l - ll2 |. Here the boundary conditions that 
4>(r > d„) is constant and <j)' {r — d a ) — Q are imposed at 
the healing distance d a , which is determined self consis- 
tently from number conservation 



(P-Pcr) 



4> 2 {d a ) 



4:TTr 2 dr = 1. 



(8) 



The prefactor p — p a = p a > takes into account that a 
given spin a only interacts and correlates with unlike 
spins a 1 ^ a. In the dilute limit 4>(r) ~ 1 and so d a = 
(9Tr/2) 1 / 3 kp 1 a , . In the unitary limit a — > ±oo the healing 

length approaches d a = (Sn) 1 / 3 ^ 1 ^ in stead. Generally 
the healing length is of order the Fermi wavelength of the 
other component, d a kF,a' ~ 1- 

For a positive scattering length the interaction energy 
A is positive and the solution to Eq. ^ is r<j>(r) oc 
sin[fc(r — b)} with X a = h 2 k 2 /2m. Defining n a — kd a the 
boundary conditions and number conservation requires 

M 



tan k„ 



1 



1 + K a tan k 



(9) 



The resulting interaction energy reproduces the correct 
dilute limit resuit of Eq. |T]) . In the unitarity limit a — S- 
+oo, the positive energy solution reduces to Ktan/î = — 1 
with multiple solutions K-y = 2.798.., K2 = 6.121.., etc., 
and asymptotically K n = mr for integer n. In addi- 
tion there is one négative energy solution for n = 
with Ko — 1.997.. which corresponds to the BCS-BEC 
crossover when a — > — oo. Generally, n = 0,1,2,.. is the 
number of nodes in the Jastrow wave function and each 
détermines a new universal limit with universal parame- 
ters depending on the number of nodes. The phase in the 
wave function is kb = tt{ii — 1/2) whenever the unitarity 
limit of n nodes is encountered. 

It should be emphasised that for positive scattering 
lengths the wave function and thus the corrélations func- 
tion between fermions of unlike spin and bosons oc r(j) ~ 
sin(fcr — 6) has a node at b/k which is somewhere within 
the interparticle distance. It does not vanish as r -s- as 
does the wave function for a short range répulsive poten- 
tial as in hard sphère scattering where a ~ R. Therefore 
the Gutzwiller approximation may well apply for hard 
sphère gases, strongly correlated nuclear fluids and liq- 
uid hélium as discussed in (H but it does not apply to the 
répulsive unitarity limit of ultracold gases when the wave 
function has to obey the short range boundary condition 
ofEq. 0. 

It is customary to define the universal function j3 — 
Ei n t /Ekin as the ratio of the interaction Ei nt and kinetic 
energy Ekm — (3/5)Ep- In the JS model the interaction 
energy per particle is £ mj = h 2 k 2 /2m — K 2 /2md 2 and 
thus (3 = (5/3)n 2 /k 2 F d 2 in the PM phase. In the FM 
phase the spin densities differ and the ratio of the average 
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Figure 2: Universal functions calculated within JS at zéro 
température vs. répulsive interaction: the ratio of interac- 
tion and kinetic energy /3, the pressure and chemical poten- 
tial and the spin susceptibility Xo/xt, ail with respect to their 
non-interactive values. Full curves include the FM transition 
whereas the dashed have FM suppressed, i.e. remain in the 
PM phase. 



interaction to kinetic energy can be considerably lower 
than p as shown in Fig. (2). 

Because Eq. ((9j) has a string of solutions for a given 
scattering length or kpa, k and (3 are multivalued func- 
tions which we distinguish by the index n — 0,1,2,... 
referring to the number of nodes in the many-body wave 
function between any two atoms [ïïl | . /3q nas been studied 
extensively in the BCS-BEC crossover and /3i in the ré- 
pulsive crossover [ÏJ-Q- In the répulsive unitarity limit 
n = 1 the universal parameter is /3i(fcj?a — > oo) = 
5k 2 /3(37r) 2 / 3 ~ 2.93. It has recently been measured for 
a 6 Li gas in two spin states [ÏJ. The chemical poten- 
tial in the optical trap almost doubles going from the 
non-interacting to the unitarity limit. Since it scales as 
p oc yl + /3i(co) we obtain /3i(oo) ~ 3 compatible with 
JS. In the following we concentrate on répulsive interac- 
tions and use j3 = p-y . 

The interaction energy for an atom with spin a dé- 
pends on the density of unlike spins and is given by 
A = K 2 /2md 2 = (3/5)Ef<t' P(a>kF,o'), where /3 is the uni- 
versal function for répulsive interactions. We obtain the 
total energy density at zéro température by adding the 
Fermi kinetic energy and the interaction energy A CT , and 
sum over particle densities (ÏÔI . Iîi1 | 



/ = E 8 " 



F, a 



P(ak F , a .)p a +fi\0) 



including a thermal free energy fx as given above. This 
expression generalizes the standard expression for the en- 
ergy density E/V = (3/5)E F (l + (3)p to finite polar- 
ization and température. The resuit can be understood 
from dimensional arguments as /? is dimensionless and 
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Figure 3: Polarization vs. ak F at T/T F = 0, 0.1, 0.2, 0.5 
from left to right. The second order transition yields a steep 
but continuous transitions P oc V ak F — ak c . The diamond 
indicates the transition point to a pure one-component (P = 
±1) FM at zéro température. 
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Figure 4: Compressibility and polytropic index vs. répulsive 
scattering length (akp) at zéro température. Both are second 
dérivâtes of the free energy and are therefore discontinuous at 
the FM transition. 



gives the répulsive energy of particles of spin a due in- 
teractions with particles of opposite spin. Note that the 
interaction energy and its dependence on polarization is 
given in terms of one universal function /3 of one variable 
only. As shown in Fig. 2 the ratio of the interaction to 
kinetic energy is reduced by the FM transition w.r.t f3. 
Expanding Eq. (fTUl) for small polarization gives 

4- = 1 + 13 + ^— P 2 + 0(P i ), (11) 
h 9 xt 

where the isothermal spin susceptibility is 
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p - -ak F P' + 



with p' = df3/d{ak F ) and P" = d 2 p/d(ak F ) 2 . In the 
dilute limit f3 = (lO/9ir)ak F and Eqs. (p} and ([12]) 
reduce to Eqs. ([2} and ((3]) respective ly. 

The spin-susceptibility calculated within JS is shown 
in Fig. 2 at zéro température. It diverges at ak F ~ 
0.90 where the universal function is P F u — 0.53. By 
equating the energy of the unpolarized gas, ~ (1+/3) with 
that of a fully polarized (one-component) gas, ~ 2 2 / 3 , we 
find that a first order transition requires (3 = 2 2 / 3 — 1 ~ 
0.59 > (3 F m, and therefore JS predicts a second order FM 
transition as shown in Fig. 1. This transition point is in 
remarkable agreement with two récent QMC calculations 
which find ak F = 0.86 [7] and ak F = 0.89 0. The 
QMC calculations could not détermine the order of the 
transition within numerical accuracy. In the BCS-BEC 
crossover a minor discrepancy was found between JS (TlT | 
and QMC [H [Tij which partly could be attributed to 
pairing which is excluded in the JS wave function. Since 
pairing is absent for répulsive interactions the JS modcl 



is expected to match the QMC calculations better near 
the FM transition. Note that the JS wave function is 
also used as a starting point in the QMC calculations of 

Refs. dm, m. 

Minimizing the free energy of Eq. ([TT|) we obtain 
the polarization P oc y— Xojxr at the onset of FM as 
shown in Fig. 3 at low températures. Full polarization 
is reached at ak F ~ 1.1 at zéro température only. 

The spin-susceptibility is related to the spin- 
antisymmetric Landau parameter as Fq = 
(m* /m)xo/XT — 1- The effective mass m* — m is 
implicitly assumed in the JS energy of Eq. ([TU)) . It 
has recently been measured in the BCS-BEC unitarity 
limit mo/m = 1.13 ± 0.03 [Ï5[ but not for the répulsive 
crossover yet. Since j3 at the FM transition point is 
comparable to \Po\i thèse two effective masses may 
be expected to be similar. The small déviation from 
m* = m only changes the universal functions and the 
phase diagram slightly at higher température. The order 
and the position of the transition is unchanged at zéro 
température. 



III. NUMBER FLUCTUATIONS 

The density or local number fluctuations have recently 
been measured in shot noise experiments for an idéal ul- 
tracold Fermi gas and by speckle noise in the BCS-BEC 
crossover [ÎH , \lï\ . The number fluctuations are measured 
in a small subvolume of the atomic cloud with almost uni- 
form density. The fluctuations in spin and total number 
of atoms are directly related to the spin susceptibility 
and compressibility respectively. 

The local fluctuations in total number can for a large 
number of atoms be related to the isothermal com- 
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pressibility kt — p 2 (dp/dfi)v.T, by the fluctuation- 
dissipation theorem 

(ANf _ 3 T KT 

N ~ 2 Tf kq [ ' 

Here, Ko = 3/(2/aEp) is the compressibility for an idéal 
Fermi gas at zéro températures. An idéal classical gas 
has kt = l/(pfcflT) such that the number fluctua- 
tions are Poisson: (AN) 2 /N = 1. The compressibil- 
ity is related to the symmetric Landau parameter as 
Fq = (m* /rn)Ko/ kt — 1. 

The compressibility can generally be expressed in 
terms of the universal function (3 [18j at zéro tempér- 
ature 

m = l + P+^ak F (3' + ±-(ak F ) 2 (3", (14) 
kt 5 10 

in the PM phase. Once the FM transition sets in, the 
ground state energy of Eq. (fTU)) is lowered due to fi- 
nit e polarization and the inverse compressibility drops 
as shown in Fig. 4 for JS. It is discontinuous when the 
second order FM transition sets in because it is a sec- 
ond derivative of the free energy which is softened by 
the spin-susceptibility term in Eq. (jlip . In the pure 
one-component FM phase the compressibility is that of 
an idéal one-component gas kq/kt — 2 2 / 3 . The peculiar 
and discontinuous behaviour of the compressibility at the 
FM transition is directly reflected in the fluctuations in 
total number according to Eq. (fl~5|) . 

If the FM transition was first order the compressibil- 
ity diverges at the phase transition, Le., k^/kt vanishes 
in part of the density région where < P < 1 (see 
Figs. 3+4). Consequently, the number fluctuation also 
diverges according to Eq. (fTB")) reflecting the density dis- 
continuity at a first order transition. 

The fluctuation-dissipation theorem also relates the 
thermal spin fluctuations to the spin susceptibility 

A(Nf — N±) 2 _ 3 T X T , . 

N 2 Tf Xo [ ' 

At the FM instability the spin-susceptibility and there- 
fore also the spin fluctuations diverge reflecting that 
phase séparation occurs between domains of polarization 
±P. Such domains were, however, not observed in the 
experiments of [1] within the spatial resolution of the ex- 
periment. 



IV. TRAP RADII AND KINETIC ENERGIES 

In experiments the atoms are confined in harmonie 
traps. For a sufficiently large number N = N a of 
particles confined in a (shallow) trap the System size R a 
is so long that density variations and the extent of possi- 
ble phase transition interfaces can be ignored and one can 



apply the local density approximation. The total chemi- 
cal potential is given by the sum of the harmonie trap po- 
tential and the local chemical potential p, a = (df / dp a )v,T 

Mo-( r ) + -jmtu 2 r 2 = ^m^R 2 . , (16) 

which must be constant over the lattice for ail compo- 
nents a — 1,2, ,.,v. It can therefore be set to its value 
at its edge R a , which gives the r.h.s. in Eq. (fTïï)) . The 
équation of state détermines the chemical potentials p a 
in terms of the universal function of Eq. ([TU)) . 

In a two-component spin-balanced System the chemi- 
cal potential and radii of the two components are equal 
(denoted p, and R in the following). In the FM phase 
their densities p(l ± P) differ but thèse FM spin do- 
mains coexist. Using the JS EoS of Eq. (fTÏÏ)) to cal- 
culate the chemical potential we can find the density dis- 
tribution from chemical equilibrium Eq. (| lé>[) including 
phase transitions and calculate cloud radii R, the root 
mean square RAIS = y/ (r 2 ) and kinetic energy E^ n — 
(k F /2m) averaged over ail particles in the trap. Thèse 
are shown in Fig. 5 normalized to their values trapped 
non-interacting ultracold atoms, R = (24iV) 1 / 6 oo , 
RM S Q = y/3/8Ro and E° m = {3/8)E° F respectively. 
Here, E° F = {hk° F ) 2 /2m and k° F = (24N) 1 ^/a are the 
Fermi energy and wave number in the centre of the trap 
for non-interacting atoms and ao = y/h/muj is the os- 
cillator length. Répulsive interactions reduce the central 
density and Fermi energy as can be seen from kp/k F 
shown in Fig. 5. As a conséquence the radii increase 
except for the RMS radius above the FM transition (see 
Fig. 5). It decreases because atoms are redistributed 
from the PM phase near the surface to the FM phase 
in the centre. The kinetic energy of the atoms has the 
opposite behaviour because repulsion increases the inter- 
action energy in the PM phase at the cost of the kinetic 
energy. 

In the récent experiment of Réf. [l| a transition is ob- 
served around ak F ~ 2.2 at températures T/E F = 0.12 
(and ak F ~ 4.2 at T/E F = 0.22). Tins transition point 
is a factor of ~ 2 larger than the FM transition point 
calculated in ail models [U-Q as well as JS. Rescaling 
ak F by a factor 2 we find very good quantitative and 
qualitative agreement with the data of llj as was found 
in the second order calculation of Réf. [6]). The distinct 
transitions in the radii, kinetic énergies and atomic losses 
are well reproduced. 

However, the FM domains were not observed in Réf. 
[ï| within the spatial resolution of the experiment. Di- 
rect observation should be possible in unbalanced spin 
Systems where macroscopic FM domain sizes can be re- 
alized. As the repulsion is increased the RMS radii of the 
minority spins increases faster than that of the majority 
spins. When the FM transition occurs the System favours 
a core with predominantly majority spins surrounded by 
a mantle of both spins in a PM phase. With increasing 
spin imbalance the majority spin purity of the FM core 
increases, i.e. the domains are effectively separated on 
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For a very elongated or cigar-shaped trap (prolate 
in nuclear terminology) , A <C 1, used in most experi- 
ments [ï^, [2(1, the collective breathing modes separate 
into a low frequency axial mode with oscillation fre- 
quency uj ax — y/3 — (7 + 1) _1 W3 and a radial mode with 

Urad = V 2 (7 + 1)^0 [13. 

The spin dipole mode is more complicated because it 
is sensitive to the spin susceptibility which diverges at 
the FM transition point. The EoS is far from polytropic 
and the délicate calculation of spin dipole modes with 
diverging spin susceptibility is beyond the scope of this 
work. The spin dipole mode is estimated within a sum 
rule approach in Réf. [22^ . 



VI. MULTICOMPONENT SYSTEMS 



Figure 5: Radius of the trapped cloud, RMS radius, kinetic 
energy and central Fermi wavenumber ail at zéro température 
and relative to their non-interaction values vs. répulsive scat- 
tering length (ak%). Dashed curved shows the RMS radius 
for a PM phase where the FM transition is inhibited. Due to 
repulsion the central density and thus k F is lower. 



a large scale. The amount of séparation and change in 
radii will dépend on the overall spin imbalance. 

Three component Systems with more than one Fesh- 
bach résonance as in 6 Li are also more complicated. For 
example, when the Feshbach magnetic field is such that 
two résonances an and are large but 023 small, the 
atoms will separate between a FM phase of 1 and a mixed 
FM phase of 2+3 with différent densities. 



V. COLLECTIVE MODES 

Collective modes have been studied intensively in the 
BCS-BEC crossover where they reveal important infor- 
mation of the équation of state (EoS) and détermine 
fixait p)- When the EoS can be approximated by a sim- 
ple polytrope P oc p 7+1 the collective eigen-frequencies 
can be calculated analytically [Ï8L [2Ï1 | in terms of the 
polytropic index 7. Even when the EoS is not a perfect 
polytrope the collective modes in the BCS-BEC crossover 
could be described well using the effective polytropic in- 
dex at densities near the centre of the trap given by the 
logarithmic derivative [ïj|] 



p dP 
= P~dp~ 



(1 + p) + \ak F p + Uak F ) 2 /3 n 



1 + P + ak F /3'/2 



(17) 



We therefore calculate 7 within JS for répulsive interac- 
tions as shown in Fig. 4. Like the compressibility it has 
a discontinuity at the FM phase transition because it is a 
second derivative of the free energy with a second order 
transition. In both the dilute limit and pure FM phase 
the gas is idéal with polytropic index 7 = 2/3. 



Interesting information on the order of the FM tran- 
sition can be obtained by generalizing the above results 
to Fermi gases with more that two spin states such as 
6 Li with v = 3 hyperfme states [Hj], 137 Yb with six nu- 
clear spin states |24j, and heteronuclear mixtures of 40 K 
and 6 Li (2f|. The interactions and phases can be very 
complicated when the Feshbach résonances between var- 
ious components differ as for 6 Li. In the following we 
restrict ourselves to multi-components with the same rel- 
ative scattering length a. 

In the dilute case the condition for a first order phase 
transition in a v component System can be found from 
the energy density of Eq. ([T]). The preferred transition 
is directly from v = 1 to a domains of one-components 
System v = 1 which occurs when 
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At zéro température this condition is akp ~ 
1.66,1.53,1.43,1.36, etc. for v = 2,3,4,5,6,.. respec- 
tively. The condition for a second order transition is 
found by expanding the dilute multi-component free en- 
ergy for small polarization. One finds the same spin 
susceptibility as in the two-component case, Eq. ((3]), 
and therefore the putative the second order transition 
remains at ak F = tt/2 ~ 1.57. Comparing numbers we 
conclude that at zéro température the second order tran- 
sition occurs for v = 2 only in the dilute case whereas for 
v > 3 the FM transition is first order and given by Eq. 
(jT5J) . At hnite températures the second order transition 
of Eq. ([3]) match the first order of Eq. ([TH)) at a tempér- 
ature which détermines the tri-critical point (akp,T) in 
the phase diagram for v > 3 as shown in Fig. 6. 

The free energy of the JS model, Eq. (jTDJ), also applies 
to multi-component Systems. The condition for a first 
order FM transition to coexisting fully polarized (one- 
component) FM domains is 
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Figure 6: Phase diagrams for multi-component (y = 
2, 3, 4, 5, 6 from left to right) Fermi gases with répulsive in- 
teractions. Full (dashed) curves indicate first (second) order 
PM to FM transitions within JS and dilute approximations. 
Circles indicate the tri-critical points where the transition 
changes from first to second order at higher températures. 



At zéro température the FM transition occurs for /3 = 
0.59,0.54,0.51,0.48,.. at akp = 0.96,0.91,0.87,0.84,.. 
for v = 2,3,4,5,.. respectively. As in the dilute case 
the spin-susceptibility is unchanged, Eq. (fl~2|). in the JS 
model and the putative second order transition remains 
when f} = 0.53 at akp = 0.90. Thus the FM transition at 
is at zéro température marginally second order for v < 3 
but first order for v > 4. Again the tri-critical points 
(akp,T) are determined by the matching condition for 
the first Eq. (fl"9)) and second Eq. (fl"2j) order transitions 
and are shown in Fig. 6. 

Generally the différence between first and second or- 
der FM transition is small which may explain why QMC 
could not détermine the order within numerical accuracy 
0, [H ■ First order transitions to partially polarized FM 
does not occur for two-component Systems but may be 
possible in multi-component Systems. 

The marginal first vs. second order FM transition 



for v = 3 is analogous to the marginal stability in the 
unitary limit of the BCS-BEC crossover [Tï| . Here it is 
known that two-component Systems are stable but four- 
component Systems are unstable as in nuclear matter. 



VII. SUMMARY AND OUTLOOK 

By extending the Jastrow-Slater approximation to fi- 
nite polarization and température we have calculated a 
number of thermodynamic functions and observables for 
cold Fermi atoms with répulsive interactions. In par- 
ticular we found a second order FM phase transition at 
akp ~ 0.90 at zéro température in close agreement with 
QMC. The compressibility and spin susceptibility were 
calculated and the resulting observables like the fluctua- 
tions in total number and spin as well as collective modes 
are discontinuous at the transition point. Thèse can be 
distinguished from a first order transition where, e.g., the 
compressibility diverges. 

For trapped gases the radii and kinetic énergies also 
have characteristic behaviour as function of répulsive in- 
teraction strength when the FM transition occurs in the 
centre. If the interaction strength is reduced by a factor 
~ 2 the radii and kinetic énergies of JS and Réf. [&] agrée 
qualitatively and quantitatively with experiments [ï[ . In 
order to observe the FM domains we suggest to start out 
with a spin-imbalanced System of two-component Fermi 
atoms and tune the magnetic field towards the Feshbach 
résonance from the repuise side where the FM transition 
sets in. As resuit the core will be a large domain of the 
majority spin only which exceeds the expérimental do- 
main size resolution. 

It would be interesting to study multi-component Sys- 
tems such as the three component 6 Li System near Fesh- 
bach résonances where bulk séparation between the spin 
component domains is predicted to take place. Multi- 
component Systems with the same interactions (and scat- 
tcring lengths) between states display interesting phase 
diagrams with first to second order tri-critical points 
when the number of components exceeds two in the dilute 
case and three in the JS model. 



Références 

[1] G-B. Jo et al., Science 32 5, 1521-1524 (2009); see also 

comment larXiv:0910.3419l 
[2] T. Bour del et al., Phy s. Rev. Lett. 91, 020402 

(2003); |cond-mat/0403091| J. Cubizolles et al., 

|cond-mat /0308018 
[3] S. Gupta et al., Science 300, 47 (2003). 
[4] E. Stoner, Phil. Mag. 15, 1018 (1933). 
[5] R. A. Duine & A. H. McDonald, Phys. Rev. Lett. 95, 

230403 (2005). 

[6] G. J. Conduit, B. D. Simons, Phys. Rev. Lett. 103, 



200403 (2009). 

[7] S. Pilati, G. Bertaina, S. Giorgini, M. Troyer, 

arXiv/1004.1169 
[8] S.-Y. Chang, M. Randeria, N. Trivedi, 

arXiv.org/1004.2680 
[9] H. Zhai, Phys Rev. A 80, 051605(R) (2009). 
[10] V. R. Pandharipande, Nucl. Phys. A 174, 641 (1971) and 
178, 123 (1971); Pandharipande, V. R., and Bethe, H. 
A., Phys. Rev. C 7, 1312 (1973); 
[11] H. Heiselberg, Phys. Rev. A 63, 043606 (2001); J. Phys. 

B: At. Mol. Opt. Phys. 37, 1 (2004). 
[12] S. Cowell, H. Heiselberg, I. E. Mazets, J. Morales, V. R. 
Pandharipande, and C. J. Pethick, Phys. Rev. Lett. 88, 



s 



210403 (2002). J. Carlson, H. Heiselberg, V. R. Pand- 
haripande, Phys. Rev. C 63, 017603 (2001). 

[13] G.E. Astrakharchik, J. Boronat, J. Casulleras, S. 
Giorgini, Phys. Rev. Lett. 93, 200 (2004). 

[14] J. Carlson, S-Y. Chang, V. R. Pandharipande, K. E. 
Schmidt, Phys. Rev. Lett. 91, 50401 (2003); S-Y. Chang 
et al., Phys. Rev. A 70, 043602, 2004; Nucl. Phys. A 746, 
215, 2004. 

[15] S. Nascimbene, N. Navon, K.J. Jiang, F. Cevy, C. Sa- 
lomon, Nature 463 (2010) 1057 

[16] C. Sanner, E. J. Su, A. Keshet, R. Gommers, Y. Shin, 
W. Hu ang, W. Ketterle , Phys. Rev. Lett. 105, 040402 
(2010); larXiv:1010.Ï874l 

[17] T. Mûller, B. Zimmermann, J. Meineke, J. Brantut, T. 
Esslinger, H. Moritz, larXiv:1005.Ô3Ô2l 

[18] H. Heiselberg, Phys. Rev. Lett. 93, 040402 (2004). 

[19] K. M. O'Hara, S. L. Hemmer, M. E. Gehm, S. R. 



Granade, J. E. Thomas, Science 298 (2002) 2179. L. Luo, 

B. Clancy, J. Joseph, J. Kinast, J. E. Thomas Phys. Rev. 

Lett. 98, 080402 (2007). 
[20] M. Bartenstein, A. Altmeyer, S. Riedl, S. Jochim, C. 

Chin, J. Hecker Denschlag, R. Grimm, Phys. Rev. Lett. 

92, 203201 (2004) 
[21] M. Cozzini, S. Stringari, Phys. Rev. A 67, 041602 (R) 

(2003). 

[22] A. Recati and S. Stringari. [ârXiv:1007.45Ô4l 

[23] S. Jochim et al, Phys. Rev. Lett. 91, 240402 (2003); M. 

Bartenstein et al., ibid. 92, 120401 (2004). T. B. Otten- 

stein, T. Lompe, M. Kohnen, A. N. Wenz, S. Jochim, 

Phys. Rev. Lett 101, 203202 (2008). 
[24] M. Kitagawa et al., Phys. Rev. A 77, 01271 9 (2008). 
[25] F. M. Spiegelhalder et al., larXiv:0908.ÎTÔTl E. Wille et 

al., Phys. Rev. Lett. 100, 053201 (2008) 



